Hall field induced magnetoresistance oscillations of a two-dimensional electron system 
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We develop a model of the nonlinear response to a DC electrical current of a two dimensional 
electron system(2DES) placed on a magnetic field. Based on the exact solution of the Schrodinger 
equation in arbitrarily strong electric and magnetic fields, and separating the relative and guiding 
center coordinates, a Kubo-like formula for the current is worked out as a response to the impurity 
scattering. Self-consistent expressions determine the longitudinal and Hall components of the electric 
field in terms of the DC current. The differential resistivity displays strong Hall field-induced 
oscillations, in agreement with the main features of the phenomenon observed in recent experiments. 
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PACS numbers: 73.43.-f,73.43.Qt,73.40.-c,73.50.Fq 



I. INTRODUCTION 



Recently, non-equilibrium magnetotransport in high mobility two-dimensional electron systems (2DES) has ac- 
quired great experimental and theoretical interest. Microwave-induced resistance oscillations (MIRO) attaining zero 
. resistance states were discovered a few years agoii^i^ in 2DES subjected to microwave irradiation and moderate 
magnetic fields, corresponding to very high Landau levels (LLs). The phototoresistance exhibits strong oscillations 
periodic in e = lu/lOc, where lu and Wc are the microwave and cyclotron frequencies. Our current understanding 
'"O ■ of this phenomenon rest upon models that predict the existence of negative-resistance states yielding an instability 
C [ that rapidly drive the system into a zero resistance stated. Two distinct mechanism for the generation of negative 
^ ' resistance states are known, one is based in the microwave induced impurity scattering^iii^i^ii^, while the second 
is linked to the modifications of the distribution produced by the microwaves and inelastic processe o^^'^^i^^ . Similar 
to MIRO are magnetoresistance oscillations induced by the combined effects of microwave irradiation and periodic 
potential modulatio n^^'^^i^^ . 

More recently. Hall field-induced resistance oscillations (HIRO) have been observed in high mobility samples in 
response to a strong DC electric curren1>ilii^. The HIRO oscillations are periodic in inverse magnetic field, with 
the resistance maxima appearing at integer values of the dimensionless parameter e = ljh/lOc- The Hall frequency 
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CS| ' ujh = li'^T^/ne) Jxl&i is associated to the energy hwn = ^R(.E?r; where the the classical Hall electric field is 
f — ■ given by = BJ^/eUe, Rc is the cyclotron radius of the electrons at the Fermi level, and 7 ^ 1.63 to 2.18. These 
\ results were confirmed in the recent experiment by Zhang etalt^, with a determination of the parameter 7 ~ 1.9. 
Additionally in this work another notable nonlinear effect was found: in the regime of separated LLs a relative 
weak DC current induces a dramatic reduction of the resistivity. Although MIRO and HIRO are basically different 
phenomena, both rely on the commensurability of the cyclotron frequency with a characteristic parameter, oj and loh 
respectively. The study of HIRO permits to analyze resistivity as a function of e when both Wc and ujh are varied 
over a wide range of frequencies. Instead the MIRO studies are carried out performing Wc-sweeps at fixed w, because 
' the experimental difficulty in implementing w-sweeps. Other interesting examples of nonlinear magnetotransport 
experiments combine DC and AC excitations^- '-P. 

Some theoretical studies of the nonlinear transport properties under DC exitations have recently appeared. In the 
work of Vavilov etal^ the nonlinear magnetotransport effects are related to changes of the electron distribution 
function induced by the DC electric field. On the other hand, the work of Lei2^ the impurity scattering processes are 
analyzed utilizing a semiclassical method in which the evolution of the CM velocity operator and relative electron 
energy are obtained from the Heisenberg operator equations, and the CM electron coordinate and velocity are treated 
classically. The aim of this work is to develop a model to describe the nonlinear response to a DC electrical current of 
a 2DES placed on a magnetic field. Our model is based on the exact solution of the Shrodinger in arbitrarily strong 
magnetic and electric fields. It is convenient to separate the electron coordinate in its relative R and guiding center X 
coordinates. Although the relative and guiding center coordinates commute; the and Ry relative coordinates satisfy 
a non-commutative algebra (and similarly for the X and Y projections of the relative coordinates). These properties 
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are exploited to work out a Kubo-like formula for the conductivity as a response to randomly distributed impurities. 
Self-consistent expressions determine the longitudinal and Hall components of the electric field in terms of the DC 
current density J^- The Hall electric field Ey is given by the dominant classical Hall contribution Ey = BJ^/erie,, 
plus an small quantum correction AEy. The differential resistance displays a strong oscillatory behavior resulting 
from current assisted electron scattering producing both intra-LL and inter-LL transitions. The results of the present 
work taken together with those of Lei^^ prove that the properties of HIRO oscillations in the overlapping LLs region 
(e > 1), as well as the strong resistance reduction in the separated LL regime can be both explained by a DC-current 
induced impurity scattering mechanism. 



II. MODEL 



The goal of this work is to study the nonlinear response to a DC electrical current density (Jx) of a 2DES placed in a 
magnetic field. The electrons are subjected both to a magnetic B = (0, 0, B) and an in-plane electric E = {Ex,Ey, 0) 
fields. Additionally the effects of impurity scattering potential V are included; hence the dynamics is governed by the 
Hamiltonian 



H^H[B,E}+V, (1) 

where 

H[B,E} = 7^J^^ + eE-x, (2) 

2m* 

includes the interaction with the magnetic and electric fields, the covariant momentum is 11 = —iKV + eA. The 
impurity potential is decomposed in terms of its Fourier components 

JVx I- ,2 

^ = E / tAi^ • {r ~ r,)] , (3) 

where is the position of the zth impurity and Ni is the number of impurities. The explicit form of V (q) depends 
on the nature of the scatterers. For short-range neutral impuritiesifl: V (q) — 2T:h^a/m* , where a is the scattering 
length, and the impurity density is related to the electron mobility according to the relation a^riimp — e/ {Air'^h^). 

Instead, in the case of a 2D screened Coulomb potential: V (q) — '^^^T'' e~'^'^ / {qxp + q), where d is the thickness of 
the doped layer and qTF = e^m* / (27reoebfi.^) ; in this case the relation of the impurity density to the electron mobility 
is approximated as Uimp = 8e{kpd)^ /{ttHij,). 

A planar electron performs cyclotron and drifting motion in magnetic and electric fields. It is then convenient 
to decomposed the electron coordinate r into the guiding center X = {X,Y), and the the relative coordinate R = 
{Rx, Ry), i-e.: r = X + R, where R = {—Hy/eB , Hx/eB). The commutation relations are 



[R.,Ry] = ^, [X,Y]=-tll, [X,,R,]=0, (4) 

with = h/eB. With this decomposition the X and Y coordinates become noncommutative (similarly for Rx and 
Ry); however the guiding center and the relative coordinates are independent variables. 

For an arbitrary orientation of the electric field E = E {cos 9, sin 9,0), the Hamiltonian H^g ^} can be exactly 
diagonalized if the vector potential is selected in the gauge: A = B {y siv? 9 , — x cos^ 9, O) . The spectrum and 
eigenstates are given by 

^^.k = , \ ^, exp{-zA:ry} exp{-(e - Cc) ^2^1} [(^ " Cc)/?!] , (5) 
V2^/i!V7r/s 

here ^ = xcos9 + ysm9, and rj = —xsm9 + ycos9 are the longitudinal and transverse coordinates (with respect to 

the electric field); = -j^ Isk, and represent the Hermite polynomials. The index k is the eigenvalue of the 

transverse —X sin 9 -\-Y cos 9 center of guide coordinate. 
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We now turn to the calculation of the current density. The velocity of the center guide coordinates are obtained 
from the Heisenberg operator equations using the total Hamiltonian in Eq. ([T]): 



Y 



- H Y] ^ + ^ 



(6) 



The current density is computed from the impurity and thermal average (J) — e(Tr [p{t) Y]) of the center of guide 
velocity V — [X, F), weighted with the density matrix that satisfies the von Neumann equation 



d 

i-g-^P=[Ho + V,p] 



(7) 



We shall now derive the Kubo-Greenwood formula for the current within the framework of the linear response 
theory with respect to the impurity potential. The electric field effects are exactly taken into account through the 
wave function solution given in Eq. The Hamiltonian is split into a unperturbed part Hs^b.e} and into the 

perturbation V{r) exp (--(5|t|); notice that we added to the impurity potential a term exp (— (5|f|), with 5 representing 
the rate at which the perturbation is turned on and off. The density is similarly decomposed as p = po + ^Pi where 
the unperturbed density matrix takes the form po = f{£a)\oL){oi\; where / is the Fermi distribution function. 
Following the usual procedure of the linear response formalism, the matrix elements of Ap, in the base given by states 
in Eq. (O, are worked out as 



(p, k\ Ap k') = k\ V \v, k') - U.k) 



(8) 



Combining the previous equations, and assuming randomly distributed impurities, the average current density is 
explicitly computed; as suggested by Eq. ^ the current is spplited in drifting and impurity scattering contributions: 



Ji{E) = —^Ej 



here the function F is given by 



j{imp) 



\V (q)r V \D^,, [/ [E^ +Lo,)-f {E,)] ImG {E^ +uo,- E,) 



IJ.,u 



and is evaluated at arbitrary values of the electric field through the argument dependence according to 



The matrix elements D^^i, are given by 



Dt,,^{q) = exp — 



UJq — el% {qyEx — qxEy 




~|2 



p>v, 



(9) 



(10) 



(11) 



(12) 



with q = (qx — iqy)/V^ and (\q\^^ denotes the associated Laguerre polynomial. 

A formal procedure to obtain the Green function requires a self-consistent calculation using the Dyson equation 
for the self-energy with the magnetic, impurity, phonon, and other scattering effects included. A detailed calculation 
of Im G incorporating all these elements is beyond the scope of the present work. Instead we choose a gaussian-type 
expression for the the density of states. This expression can be justified within a self-consistent Born calculation that 
incorporates the magnetic field and disorder effects^^i^iiS^iS^, hence the density of states for the /i-Landau level is 
represented as 



Im G (wq + E^- E^j) = \/ exp 



Eu - E^Y 



2r2 



^2 ^ m'^c 



Ttr 



(13) 
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The parameter [3 depends on the difference of the transport scattering time Ttr obtained from the mobihty and the 
single particle scattering time Tg. In the case of short-range scatterers Ttr — Tg and /3 — 1. In the case of long-range 
screened potential, Pi, depends on the filhng factor v, e.g. (3^=50 ~ 13.5.^" 

The current in Eq. ^ applies in general in the non-linear transport regimen, for an arbitrary strength of the DC 
electric field. The linear response regime is recovered if the function F{ujq) is expanded to first order in E (the zero 
order term cancels because of the angular integration) to yield: Jx — Oxx Ex and Jy — ayy Ey, where 

^..=<^yy = elsj j^,y [^J ^^^^ (14) 

and ffxy = (Tyx = en/ B. 

Let us now consider the nonlinear transport regime. In a typical experimental configuration the electric field is not 
explicitly controlled, instead the longitudinal current Jx is fixed to a constant value, while the transverse Hall current 
Jy cancels. Consequently Eqs. ([9]) lead to the conditions 

JX= ^Ey + j(/"^^\Ex,Ey), 

= -^Ex + J^^"^^HEx,Ey). (15) 

They represents two implicit equation for the unknown Ex and Ey] the equations can be solved following a self- 
consistent iteration. However, it is easily verified that for the conditions that apply in recent experiments, the 
solution of the previous equations simplify assuming that the following conditions Ex <C Ey and J^*"'^-' ^ eUeEy/B 
are simultaneously satisfied. It is then possible to explicitly solve for Ex and Ey 

Ex = —Jr^Ex, Ey) « —J"^P{Ex =0,Ey = BJx/en,) . (17) 
erie eue 

In order to derive the previous results we used the fact that J™'^p{Ex — Q,Ey) cancels because of the angular 
integration. Eq. (fT6|) shows that the leading contribution to the Hall electric field is given by the classical result 
= BJx/eUe, however there is a correction AEy given by the second term of the RHS. of Eq. ITBl) . Utilizing Eqs. 
Q the correction to the Hall electric field can be worked out as 



hEx f cPq 2 f dF 
eue J C 

where Ex is determined by Eq. p7p and uj* is given by 



QxJx 



eUe 



(19) 



Figure [T] shows the electric fields: Ey = BJx/eUe, AEy, and Ex as a function of the magnetic field for a fixed 
value of the longitudinal current. Notice that the conditions AEy ^ Ex ^ Ey'' hold in general, validating the 
assumed approximation. In fact it is verified that the approximated solutions in Eqs. pop and (|17p coincide with the 
self-consistent solutions that are obtained from Eqs. with a better that 1% precision. 

Collecting the previous results, it follows that the Hall resistivity is well approximated by the expression 



p^y^En^^^lE^f^) (20) 

Jx en^Jx \ 9Ex J {^E^=0,Ey=BJ^/en^) 



Whereas the expression for the non-linear longitudinal resistivity is given as 

Ex ^ B r J^q 

Jx eUeJx J (Stt 



= ^ - -z:-^ I TTd^^-FH) ' (21) 
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with u!* given in Eq. (fTQ]) . The differential resistivity is calculated as r^x = d{JxPxx) /dJx yielding 



ieue 



Eqs. ([9]) for the nonlinear current together with the definitions in Eqs. (|f Ollf 3[) constitute the central result of 
the paper. They apply in general in the nonlinear regime in which both the longitudinal and Hall electric fields are 
arbitrarily strong. However, for the conditions that apply in the experiments of current interest, it is reasonable to 
consider the Ex weak limit. Then, the Hall field is accurately approximated by the classical result Ey = BJx/eUe, 
whereas pxx and the differential resistivity are explicitly computed from Eq. (j2ip and Eq. (|22p respectively. 

1/2 

In the work of Zhang etal^^ the Hall frequency is defined ujh = 'yJx (27r/e^ne) . Here we assume that 7 = 2, 
this can be justified if we observe that the integral in Eq. (j2ip is evaluated in terms of the variable uj* = qxJx/en,, 

and it is dominated by contributions of exchanged momentum in the region qx ~ 2kp. Recalling that kp = \/2Trn^, 

1/2 

it yields ujh = Jx (STr/e^rie) . The dimensionless control parameter is then given by the ratio of the Hall to the 
cyclotron frequencies 

UJH 2eEyRc ^ vp 



UJr huJr UJr 



Rc = —, (23) 



and it can be interpreted as the ratio of the work of the electric Hall field associated with the displacement of the 
guiding center of the cyclotron trajectory by 2Rc to the Landau energy hujc- 



III. RESULTS 



Results are presented for a 2DES GaAs sample, with parameters corresponding to those reported in recent experi- 
ment of Zhang etal^: m* = 0.0635TOe , electron mobility p, — 1.2 x 10^ cm? /V s and density rig = 3.7 x lO^^cm"^, 
and lattice temperature T = 1.5K. For the impurities we consider an average of short- and long-range scatterers, 
selecting the parameter f3 that appear in the Eq. (fT3|) as /3 = 5. 

Fig. O (a) shows the differential resistance rxx as a function of the magnetic field B for a fixed current density 
Jx — Q.8A/m {ujh/27t « 65GHz). We observe clear differential magnetoresistance oscillations mounted on an offset 
of 1.4iri determined by the Drude contribution. At the top of this figure the values of e = ujh/uJc are displayed, 
suggesting an oscillations period Ae ~ 1. To confirm these observations, in Fig. [l](b) the Hal-field induced correction 
^Txx = rxx — rxx{Jx = 0) is plotted as a function of e. Magnetoresistance oscillations are clearly observed up to the 
seventh order. The first peak appears at e '--^ 0.95, for higher e oscillations the maxima occur at e « j, with j an 
integer; while the minima are very close to e ~ j + 1/2. These results are very similar to the experimental findings of 
Zhang etal:^, although the localization of maxima (minima) close to integer (half integer) are here obtained when 
7 = 2, whereas in that work they correspond to the selection 7 = 1.9. The amplitude of the differential resistance 
oscillations display a rapid decay as the magnetic field decreases. This decay can be parametrised by the Dingle 
factor S = exp {—n/uJcTs), this allow us to get an estimate of the single particle scattering time « I5ps, that 
compares with the transport scattering time as Ts/rtr ~ 4.8 ; in good agreement with the selected value /3 = 5 . 

To further support the previous results, we present in Figure [3] (a) a plot of Vxx as a function of l/B for various 
values of the longitudinal current density Jx- As expected each curve shows a different period as a function of l/B, 
though the decay of the oscillation amplitudes are well described by the same Dingle factor d = exp {—tt/uJcTs), with 
Ts « 15 ps. Notice that the all the curves share a maximum a,t B = 1/2.5KG given that OAA/m, 0.6 A/m and 
0.8 A/m are multiples of 0.2 A/m. Notwithstanding, when /S.rxx is plotted as functions of e in Fig. [3] (b), the positions 
of the maxima (minima) of the various curves draw close to integer (half integer) values of e. 

The oscillation waveform can be more clearly appreciated in graphs in which e varies performing current-sweepes 
at fixed B. In [4] (a) and (b) plots of /S.rxx (e) are displayed for B = 416 G and B = 732 G. We again observe that 
the first maxima appears at e < 1, but as e increases, the position of the maxima (minima) are localized very close to 
integer (semi-integer) values of e. The oscillation amplitude is almost constant, because the Dingle factor now remains 
constant as e varies. 

We finally focus our attention to the regime of separated LLs, e < 1. Fig. [5] (a) displays the differential resistivity 
as a function of the magnetic field for both zero DC bias and a small current Jx = 0.25 A/m. Notice that for the 
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selected small current [J^ = 0.25 ^/m) only one HIRO peak is resolved, however in the region above above this peak 
{ujh < LJc) the resistivity is strongly suppressed. This suppression is in very good agreement with the experimental 
results of Zhang eioL— , and resembles the observed suppression of resistance observed for MIRO. The suppression is 
also observed in[5](b) where Jx— sweeps are performed in order to plot Ar^xi^) at fixed values of B; several values 
of B are selected, from top to bottom B = 0.75, 1.0, 1.25, 1.5, 1.75, 2.0 fcG. It is observed that the minima becomes 
deeper as B is increased, but in all the cases the the maxima lies very near e = 1. The values of the width F (Eq. 
(fT5)) ) are marked in the plot by the vertical dotted {B = 0.75 kG) and dashed {B = 2.0 fcG) lines. Considering the 
expression for the density of states in Eq. p3p . we can interpret the decaying of pxx m the region uh < T as a result 
of intra-Landau transitions; whereas in the region 2 F < luh the resistivity increases as a result of the inter-Landau 
transitions. A similar behavior is observed when the resistivity is plotted as a function of e; however now the minima 
of 

Pxx is attained at e — 0.6 for all values of B. The minima of Pxx corresponds to the region 2 F < lujj < ^ 2 F in 
which both intra-Landau and inter-Landau transitions are both suppressed. 



IV. CONCLUSIONS 



In conclusion we presented a theory for the non-linear transport of a two 2DES placed in a magnetic field. The non 
linear response to a DC current is incorporated by the exact solution of the Schrodinger equation including the effects 
of arbitrarily strong magnetic and in-plane electric fields. By means of the non-conmuting relative and guiding center 
coordinates, a linear Kubo formula with respect to the impurity scattering is worked out. The nonlinear expression 
for the electric current Eqs. ((9] |lip constitute the central result of the paper. They apply in general in the nonlinear 
regime in which both the longitudinal and Hall electric fields are arbitrarily strong. However in the experiments of 
current interest, it is reasonable to consider the Ex weak limit in which the Hall field is accurately approximated by 
the classical result Ey = BJx/eUe, thus pxx and the differential resistivity Vxx are explicitly computed from Eq. (|2ip 
and Eq. (|22[) respectively. Our model is able to reproduce the most important features of recent experiment a^^i^^'^^ . 
In the region of separated LLs the differential resistance as a functions of e presents strong oscillations with constant 
period (e = 1); the dominant mechanism being the current-induced impurity tunneling between Landau levels. In 
the region of of separated LLs the dramatic reduction of the resistivity at relative weak field is well reproduced; the 
origin being related to the suppression of scattering within the LLs. 
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FIG. 1: (color onhne) Induced electric Hall field: classical E"^ (B) and quantum correction AEy (B), and longitudinal electric 
field Ex (B)) calculated from Eqs. (|16|l and (|17p at Jx = 0.8 A/m. The values of the other parameters are: m* = 0.0635me, 
= 1.2 X 10'' cm^/Vs, Tie = 3.7 x lO^cm-^, T = 1.5K, and (3 = 5. 



8 




FIG. 2: (color online) (a) Differential resistance r^x (B). (b) Correction to differential resistance Ar^x = r^x — rxx{Jx ~ 0) 
versus e — uih/^^c- The parameters have the same values as in Fig([l|. 
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FIG. 3: (color online) (a) Differential resistance r^x as a function 1/B and (b) Correction to the differential resistance Ar^x (e) 
for: Jx ~ OAA/m (red), Jx ~ 0.6 A/m (blue) and Jx = 0.80 A/ m (green, higher peak). 
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FIG. 4: (color online) Correction to differential resistance l\rxx (e) obtained from J^^-sweeps at fixed values oi B = 416G (a) 
and B = 712G (b). 
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FIG. 5: (color online) (a) Differential resistance Vxx (B) for Jx ~ (black) and Jx ~ 0.25A/m (red), (b) Correction to 
differential resistance Arxx (e) obtained varying Jx at fixed values of B, from B = 0.5 kG (red) to B = 1.75 kG (light blue) 
in 0.25 fcG steps (top to botton at e = 0.4). (c) Normalized longitudinal resistance as a function of uh for the same selection 
of values of the magnetic fields as in (b). The vertical line shows the value of 2T, Eq. (|13p . for B — 0.5 kG (dotted) and 
B — 1.75 kG (dashed) showing good agreement with the half- width of the Jx decay of pxx/pxx{0). (d) This figure displays 
similar plots as a function of e. 



